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ABSTRACT 

The extended Baryon Oscillation Spectroscopic Survey of the Sloan Digital Sky Survey (SDSS-IV/eBOSS) has an 
extensive quasar program that combines several selection methods. Among these, the photometric variability technique 
provides highly uniform samples, unaffected by the redshift bias of traditional optical-color selections, when z = 2.7 —3.5 
quasars cross the stellar locus or when host galaxy light affects quasar colors at z < 0.9. Here, we present the variability 
selection of quasars in eBOSS, focusing on a specihc program that led to a sample of 13,876 quasars to (jdered = 22.5 
over a 94.5 deg^ region in Stripe 82, an areal density 1.5 times higher than over the rest of the eBOSS footprint. 

We use these variability-selected data to provide a new measurement of the quasar luminosity function (QLF) in the 
redshift range 0.68 < « < 4.0. Our sample is more dense, and reaches deeper than those used in previous studies of 
the QLF, and is among the largest ones. At the faint end, our QLF extends to Mg{z = 2) = —21.80 at low redshift 
and to Mg{z = 2) = —26.20 at a ~ 4. We fit the QLF using two independent double-power-law models with ten free 
parameters each. The hrst model is a pure luminosity-function evolution (PLE) with bright-end and faint-end slopes 
allowed to be different on either side oi z = 2.2. The other is a simple PLE at 2 ; < 2.2, combined with a model that 
comprises both luminosity and density evolution (LEDE) at 2 > 2.2. Both models are constrained to be continuous at 
2 = 2.2. They present a flattening of the bright-end slope at large redshift. The LEDE model indicates a reduction of 
the break density with increasing redshift, but the evolution of the break magnitude depends on the parameterization. 
The models are in excellent accord, predicting quasar counts that agree within 0.3% (resp., 1.1%) to g < 22.5 (resp., 
g < 23). The models are also in good agreement over the entire redshift range with models from previous studies. 

Key words. Quasars: general, large-scale structure of Universe, surveys 


1. Introduction 

Quasars have become a key ingredient in our understand¬ 
ing of cosmology and galaxy evolution. Being among the 
most luminous extragalactic sources, they have become a 
mainstay of cosmologic al surveys such as t he 2dF Quasar 
Redshift Survey (2QZ; Groom et al. 2001 1 and the Sloan 
Digital sky Survey (SDSS; York et ai. 20001, where they 
are the source of choice to study large-scale structures at 
high redshift. Quasars can be used as direct tracers of 
dark matter in the redshift range 0.9 < z < 2.1 where 
they are present at sufficiently high density, and as back¬ 
ground beacons to illuminate the intergalactic medium 
at higher redshift, where the cosmological information is 
produced by the foreground neutral-hydrogen absorption 


systems that form the Lyman-a forest. As part of the 
thi rd-generation of the S loan Digital Sky Survey (SDSS- 


III; Eisenstein et al.|2011[), the Bar yon Acoustic Oscillation 


Survey (BOSS; Dawson et al.|2013 1 measured the spectrum 


of about 300.000 quasars, 180.000 of which are at z > 2.15, 
to a limiting magnitude of 5 ~ 22. As part of SDSS-IV, the 
extended Baryon Oscillation Spectroscopic Survey (eBOSS ; 


Dawson et al. 2015 Tinker & SDSS-IV Collaboration 2015) 
IS aiming to more than quadruple the number of known 
quasars over redshifts of 0.9 < z < 2.2 to 5 ~ 22, in addition 
to targeting new quasars at z > 2.2. The next-generation 
Dark Energy Spe ctroscopic Instrume nt (DESI, previously 
named BigBOSS; Schlegel et al.|20lT I is designed to obtain 
spectra of more than two million quasars, reaching limiting 
magnitudes g ~ 23. This new challenge requires, as a first 
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step, a good knowledge of the quasar luminosity function 
(QLF) in order to determine the expected number count 
for quasars, and optimize the distribution of fibers among 
the various cosmological probes. 

In the past two decades, with the advent of large quasar 
surveys, the number of known quasars has increased by over 
a factor 20, triggering significant effort to measure the QLF 


(see Ross et al. 2013, for an overview of recent determi- 
nations). Nevertheless, the measurement of the QLF over 
2 < z < 4, where the number density of quasars starts to 
decline and their selection with traditional color-based al¬ 
gorithms is less efficient, remains challenging, especially at 
the faint end. This situation arises because the broad-band 
colors of z ~ 2.7 and z ~ 3.5 quasars ar e, respectively. 


very similar to those of A-F and K stars (|Fan 


1999 Fan 


et al. 2001 Richards et al. 2002 Ross et al. 2012[). The 


2005), 


at 


Oscillations, or BAO (Eisenstein et al. 
level in the yet unexplored 0.9 < z < 
plans to target and spectroscopically 
500,000 quasars in this redshift range, including quasars 
already confirmed with SDSS-I/II. At z > 2.1, eBOSS 


the 2% 
regime, eBOSS 
identify at least 


will complement previous studies from BOS S ( Slosar et al.^ 
2013 Busca et al. 12013 Delubac et alj201^ I and provide a' 
measurement of the BAO feature at 1.57o, using a sample 
of 75,000 quasars that had not previously been identified. 
In addition, eBOSS has conducted an extensive search for 
quasars at all redshifts in a 120 deg^ area where unique 
time-domain photometry from SDSS is available. Because 
this region allows a highly complete selection of quasars 
with minimal completeness corrections, it is ideal for QLF 
studies, and is the focus of the prese nt work. We improv e 
upon previous QLF studies^ such as Groom et al. ( |2009| ) 
Palanque-Delabrouille et al. (2013) and" Ross et al. 




in terms of the size of the quasar sample used for tJie mea¬ 
surement, in depth, and in redshift homogeneity of the tar¬ 
get selection. For the quasar luminosity function, we build 
on earlier semi-empirical models (e.g., 


Schmidt & Green 


et al.| 2004 Richards et al.| 2005 

dasinger et al. 2005 

Richards et al. 

0 

0 

CN 

Hopkins et a 

2007 

) such as Pure 


tially with look-back time, Luminosity Dependent Density 
Evolution (LDDE) and Luminosity Evolution -|- Density 
Evolution (LEDE). 


density of z <0.9 quasars is not well-characterized either 
as host galaxy light can significantly affect the colors of 
faint quasars. To circum vent these difficulties, [Palanque^ 
Delabrouille et al. (2011) developed a selection algorithm 
relying on the time variability of quasar fluxes. This tech¬ 
nique was demonstrated to increase by 20 to 30% the den¬ 
sity of identified quasars, and to effectively recover addi¬ 
tional quasars in the redshift range 2.5 < z < 3.5. It 
was applied to measure the QLF over the redshif t range 
0.68 < z < 4.0 (Palanque-Delabrouille et al. 2013), using 
a sample of quasars that was found to be 80% complete 
to g = 20, and still 50% complete at 5 = 22.5. Despite its 
limited statistics of 1877 quasars, this study yielded com¬ 
petitive results which have been used to estimate quasar 
counts for several ongoing large-area surveys. 

The QLF can only be improved with a well-controlled 
quasar sample of much larger size. In the near-term, eBOSS 
is the most ambitious survey verifying this requirement. 
Designed to measure the scale of the Baryon Acoustic 


In this paper, we present a sample of 13,876 quasars, se¬ 
lected by eBOSS over a 94.5 deg^ region with a technique 
relying upon quasar time-domain variability. For this study, 
we take advantage of spectroscopy conducted by eBOSS of 
the part of the SDSS So uthern equatorial stripe , hereafter 
referred to as Stripe 82 (Stoughton et al. 2002), where 50 
to 100 epochs of imaging are available over a time period 
of about 10 years. The variability technique used here is a 
robust, efficient and well-understood method whose com¬ 
pleteness can be readily evaluated using an independent 
control sample. With this strategy, all completeness correc¬ 
tions can be derived from the data, without requiring any 
model of quasar light curves or colors. 

The outline of the paper is as follows. In Sec. we 
present the variability programs in eBOSS. In Sec. we 
describe the imaging data, provide the details of the selec¬ 
tion of the targets for this study and present the resulting- 
spectroscopic data. In Sec. we give the raw quasar num¬ 
ber counts, explain the computation of the completeness 
corrections, and derive the completeness-corrected number 
counts. Finally, in Sec. we derive the QLF from our 
data. The present analysis refers extensively to our previous 
works on quasar variability. To simplify the presentation 
and make it easier for the reader to identify any references 
to these earlier papers, we will henceforth refer to our paper 
demonstrating the use o f time-domain variability for quasa r 
selection as Paper Var ( Palanque-Delabrouille et al.||2011 ), 
and to our paper pres enting our previous measurement o f 
the QLF as Paper LF (Palanque-Delabrouille et al.||2013|. 


2. Time-domain quasar selection with eBOSS 

Data for SDSS-IV/eBOSS is taken with the 2.5-meter Sloan 
Foundation Telescope (Gunn et al. 2006), using the same 
spectrograph and data reduction pipeline as for SDSS- 


III/BOSS (Bolton et al. 2012 Smee et al. 2013 


et al. 2013). The eBOSS surve} ^Dawson et al. 


eludes an extensive quasar program ([Myers et al. 


Dawson 


2M^ 


m- 

. A 

CORE selection will provide a homogeneous sample of at 
least 69 deg“^ 0.9 < z < 2.2 quasars, and the combination 
of several techniques will increase the sample of z > 2.1 
quasars by more than ~ 7 deg~^ compared to BOSS. The 
majority of the quasars at z > 2.1 are obtained either from 
the GORE selection, which is not strictly limited to z < 2.2 
and provides of order 6 deg~^, or from a selection based on 
quasar variability that provides another ~ 3 deg~^ quasars. 

The use of time-domain photometric measurements to 
exploit quasars’ intrinsic variability has been demonstrated 
during the course of the BOSS survey in Paper Var and 
Paper LF. In the context of eBOSS, variability selection 
of quasars is performed over 90% of the survey foot¬ 
print using tim e-domain data from the Palomar T ransient 
Factory (PTF: 


Rau et al. 2009 Law et al. 2009). Details 


on the variability selection from RTF data are available in 


Myers et al. (20151. Over Stripe 82, however, the SDSS pro- 
vides data that are both deeper and with a longer lever-arm 
in time than PTF. In this region, we therefore replace the 
PTF selection by a dedicated program (VAR_S82) based on 
a variability selection of quasars from SDSS photometry. 

In Tab. we list the three selection methods dedicated 
to quasars, with their eBOSS targeting bit names, numer¬ 
ical equivalents, and average target density over Stripe 82 
for GORE and VAR_S82, and over the eBOSS footprint 
where PTF data are used, thus outside Stripe 82, for PTF. 
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We also provide the density of quasars already identified 
spectroscopically in Stripe 82 (hereafter referred to as the 
‘known’ quasars), which is higher than over the rest of the 
eBOSS footprint because of several BOSS programs dedi¬ 
cated to quasar selection in Stripe 82. The listed density for 
CORE is after removal of the overlap with known quasars, 
and the density for both RTF and VAR_S82 are given after 
removal of the overlap with both known and CORE sam¬ 
ples. 


Bit 

Name 

Density (deg ^) 

- 

Known quasars 

80 

2io 

QSO_EBOSS_CORE 

60 

2“ 

QSO_PTF 

20 

2® 

QSO VAR S82 

50 



Table 1: eBOSS quasar targeting bits and average targeting 
densities in Stripe 82, except for PTE. The CORE density 
is that after removal of the overlap with known quasars. 
The PTE and VAR_S82 densities are given after removal of 
the overlap with both CORE and known quasars. 


All BOSS quasar targets were visually inspected to b e 
classified as star, galaxy or quasar ( Paris et al.|2012 2014). 
This procedure is no longer possible in eBOSS where the 
density of quasar targets is increased by at least a fac¬ 
tor 4 (about a factor 5 in specific regions such as Stripe 
82). Studies on the pipeline performance allowed an im¬ 
provement of the consistency between pipeline and visual- 
inspection classifications, and thus significantly reduced the 


need for visual inspection (cf. Dawson et al. 


2015 for de¬ 


tails on how spectra are selected for visual inspection). We 
identified the types of failures that could not be systemat¬ 
ically associated with any given class. The remaining fail¬ 
ures are flagged as requiring visual inspection. As a result 
of these improvements, the fraction of visually-inspected 
quasar targets, including both those that were identified as 
needing an inspection and those belonging to the few ran¬ 
dom plates that were visually inspected for evaluation of 
the pipeline performance, is of order 8% on average over the 
eBOSS footprint. For the QSO_VAR_S82 targets, however, 
this fraction rises to ~17% due to the fainter brightness on 
average of the selected objects. 

For all the objects, the eBOSS pipeline encodes the 
spectrum classification into the CLASS_AUTO flag, and 
the redshift, when relevant, into Z_AUTO. When the spec¬ 
trum was visually inspected, two additional flags are set: 
CLASS_PERSON, which encodes the classification, and 
Z_CONF_PERSO N, which eii codes the confidence on the 
redshift estimate ( Paris et'n(]|2014 ). We define a spectrum 
“uberclass” as follows. If a visuai inspection was done and 
led to a clear identification (Z_CONF_PERSON> 2), then 
the object uberclass is set equal to CLASS_PERSON (i.e.. 
Star, Quasar or Galaxy). If visual inspection did not lead to 
a clear identification, the uberclass is set to ‘Inconclusive’. 
In the absence of visual inspection, the uberclass is set equal 
to CLASS_AUTO, which can also be Inconclusive. We de¬ 
fine our sample of quasars as the set of targets with uber¬ 
class equal to Quasar. The quasar redshift is set equal to 
the visual inspection redshift if the latter is available, and 
to Z_AUTO otherwise. 



Fig. I: Redshift and magnitude distributions of the quasars 
selected over Stripe 82 with the QSO-EBOSS-CORE (lOfSl 
objects) or the QSO-VAR_S82 (16243 objects) flag. 


We illustrate in Fig. [T] the magnitude and redshift dis¬ 
tributions of the quasars selected in Stripe 82 with the 
QSO_EBOSS_CORE or the QSO_VAR_S82 flag, includ¬ 
ing the known quasars (and the overlap with CORE in 
the variability sample). A large fraction of the quasars 
are common to both selections. The QSO_VAR_S82 sam¬ 
ple, however, contains about 1.6 times more quasars than 
the QSO_EBOSS_CORE sample. The origin of this im¬ 
provement is 2-fold. Part of it is due to the fact that 
the QSO_VAR_S82 sample is selected from about 50 
epochs of photometry, instead of a single epoch for the 
QSO_EBOSS_CORE sample (which is done to ensure the 
uniformity with the rest of the eBOSS footprint). The other 
part of the improvement comes from the different selection 
techniques: at identical depth for the input photometry, us¬ 
ing for instance 50-epoch coadded images to measure object 
colors and 50 individual epochs of imaging to measure vari¬ 
ability criteria, we have shown in Paper Var that the vari¬ 
ability selection selects 30% more quasars than the CORE 
selection for the same total number of targets. 

The QSO_VAR_S82 sample significantly increases the 
completeness at all redshifts and magnitudes, and in par¬ 
ticular at the faint end {g ~ 22-22.5). In the rest of this 
paper, we focus on the targets with the QSO_VAR_S82 bit 
set. 
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3. Data and target selection 

This section provides an overview of the control sample of 
quasars used throughout this study. We present the imag¬ 
ing data from which the targets are selected, describe the 
selection algorithm and discuss the spectroscopic survey. 
The magnitudes of the sources are denoted rt, g, r, i and z 
when referring to observed magnitudes, and Udered, Sderedj 
Cdered) Mered and Zdered when referring to magnitudes cor¬ 
rected for G alactic extinction usi ng the extinctions from 
the maps of Schlegel et al. ( 1998|). The bands correspon d 
to the SPSS filters ( [Fukugita et ^|1996 Doi et al^|2010 |. 


3.1. Control sample 

The completeness corrections related to the criteria used 
to select the targets are determined using a control sam¬ 
ple of 4555 spectroscopically-confirmed quasars that were 
selected in Stripe 82 independently of any variability cri¬ 
terion. This prevents our control sample from being biased 
towards sources that exhibit a high quasar-like variability. 
Such quasars would indeed have led to over-optimistic com¬ 
pleteness estimates for our selection algorithm. 

The control sample is built from the 2dF quasar cata¬ 
log (Groom et al. 20041, the 2dF-SDSS LRG and Quasar 
Survey 2 SLAQ (G room et al. |2009 ), the SDSS-DR7 quasar 
catalog ( Schneider et al. |2010[) and BOSS observation s 
through August 2010 (Ahii et al. 2012 Paris et al. 2012). 
These catal ogs were obtained from pure color se lections (cf. 
for instance iRichards et al.l (2002 1 for DR7 and Ross et al. 

obi 


(2012) for BOSS). BOSS observations taken after Summer 
2010 on Stripe82 had contributions from a variability selec¬ 
tion and were therefore discarded from the control sample. 

Golor information comes from flux ratios, thus not sensi¬ 
tive to absolute fluxes, while variability information comes 
from variations in absolute fluxes. Furthermore, time vari¬ 
ations are seen to be synchronous in different bands, thus 
not affecting the source colors. Color and variability selec¬ 
tions are therefore complementary, as was already shown 
in Paper Var, with no obvious correlations between the 
photometric and time-domain characteristics of quasars. 
The only source of correlations could come from the im¬ 
age depth, but we take this into account by computing all 
corrections as a function of source magnitude. We thus ex¬ 
pect no measurable bias in our completeness estimates from 
the use of this control sample. 

з. 2. Imaging data and target selection 

The selection of the targets for this study relies heavily on 
the variability selection described in Paper Var where all 
the details can be found; we therefore only summarize the 
main steps here. 

The initial source list is det ermined from the c o-addition 
of single-epoch SDSS images (Annis et al. 2014[ ) in Stripe 
82, from which we take the source magnitude (in SDSS 

и, g, r, i and z bands) and morphology. As morphology 
indicator, we use a continuous variable defined as 


rridiS = rupspig) ~ ”^model(5) ) 


( 1 ) 


where mpsF(ff) and mi„odei(ff) are the magnitudes of the 
source, measured in the g band, obtained from a PSF fit 
(valid for unresolved objects) or a model fit (more ap¬ 
propriate for extended objects, where model can be a de 


Vaucouleurs or an exponential shape, for instance), re¬ 
spectively. As shown in Fig. mdiff peaks near 0 with a 
standard deviation of 0.01 for point-like sources, and ex¬ 
tends from 0.01 to values beyond 0.3 for extended sources. 
Because the emitting region of a quasar, too small to be re¬ 
solved, outshines the host galaxy by a large factor, a quasar 
generally appears as a stellar-like source. Only for the near¬ 
est quasars (redshift z < 0.9 at most) can the host galaxy 
be detectable in co-added images, making the source ap¬ 
pear extended in those cases (see Sec. 2.4 of Paper LF or a 
quantitative study of the effect). The average morphology 
indicator as a function of redshift for the quasars of the 
control sample is displayed in Fig. It is at the level of 0.2 
or more as the redshift approaches 0 , has a value of about 
0.1 at z ^ 0.5, and is below 0.05 for redshifts z > 0.6. At 
z > 2.5, the small decrease in mdiff with increasing red¬ 
shift is correlated with the increasing average magnitude of 
the quasars. The larger the photometric errors, the less the 
PSF and model fits tend to capture a difference in spatial 
extension of the source, and hence the closer on average is 
mdiff to zero. We apply a cut on morphology (upper bound 
on mdiff, cf. details below) to reject galaxies. 

To reduce the stellar contamination, we apply a loose 
color cut by requiri ng that C 3 < 1 — 0.33 x ci, where ci and 
C 3 are defined as in Fan (1999) by 


Cl = 0.95(u — 5 )-I- 0 . 31(5 — r)-f 0.11(r — i), 
C3 = —0.39(m — g) + 0 . 79(5 — + 0.47(r — i). 


( 2 ) 


This is the same criterion as was used in Paper LF, where 
we had estimated that close to 100% of z < 2.2 and 98% 
of z > 2.2 k nown quasars (i.e., quasars from the control 
sample of Sec. 3.1) passed that condition. The completeness 


of this color cut is in cluded in the selection efficiency Csei 
described in Sec. 14.11 

The main selection is based upon a criterion measuring 
the variability with time of the source. The lightcurves of 
our sources contain, on average, 52 SDSS individual epochs 
spread over 7 years. They are used to compute two sets of 
parameters that characterize the source variability: 

- the of ffio fif of the lightcurve in each of the ugriz 

filters by a constant m: ~ w)/(Ti]^, where the 

sum runs over all observations i, 

- two par ameters, an amplitud e A and a power 7 as intro¬ 
duced by Schmidt et al. ( 2010 |, that characterize the vari¬ 


ability structure function V(Atij), i.e., the change in mag¬ 
nitude Amy as a function of time lag Afy for any pair ij of 
observations: V(Atij) = |Amij| —(crf-l-cr?)^/^ = Ax(Atij)^. 
Because quasars have similar time variations in different 
bands, we reduce the uncertainty on variability parameters 
by fitting simultaneously the 5 , r and i bands (those least 
affected by noise and observational limitations) for a com¬ 
mon 7 and independent amplitudes Ag, and Ai. 

Variable objects, whether quasars or stars, are expected 
to have large y^’s, thus allowing a distinction between 
variable and non-variable targets. The structure function 
parameters A and 7 can discriminate between these two 
classes of variable objects: quasars tend to have both large 
A and large 7 , due to magnitude changes that increase with 
time, while variable stars (such as pulsating or eclipsing bi¬ 
naries) can have large A but usually 7 near 0. 

For each source, a neural network combines the five 
X^, the power 7 and the amplitudes Ag, A^ and 
to produce an estimate of quasar-like variability. The 
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mag PSF - mag Model (g) 

Fig. 2: Difference between PSF and model magnitudes, used 
as morphology indicator, for random objects in Stripe 82. 
Point-like sources in the co-added images, overlaid in blue, 
have a small magnitude difference. Quasar selection re¬ 
quires magnitude differences of at most 0.05, relaxed to 0.1 
for objects with high significance of quasar-like variability 
(to recover of low-z quasars where the host galaxy can make 
the source appear extended). 



Redshift 


Fig. 3: Morphology indicator as a function of redshift for 
quasars. The cut at 0.05 includes most quasars at z > 0.6, 
but rejects most lower-redhift quasars whose host galaxy is 
detectable in co-added images. 


(unn > 0.85), a strict morphology cut (wdiff < 0.05) is ap¬ 
plied in case of a lower variability indicator (0.50 < ^nn < 
0.81), and for values of unn intermediate between 0.81 and 
0.85, the threshold on muiff is gradually evolved between 
the two extreme values of 0.05 and 0.1. Even the tightest 
bound (0.05) fully encompasses the range of potential val¬ 
ues for point-like sources. 

We restrict the study to sources with g < 22.8 and 
5 dered < 22.5. The average g Galactic extinction over the 
observed zone of Stripe 82 is 0.12, so both limits are com¬ 
parable. With this magnitude limit, the selection described 
above leads to a sample density of 175 deg“^ targets. 
Removing targets that already have a spectroscopic identi¬ 
fication from previous observations reduces the sample to 
about 95 deg~^ targets. Further removing the overlap with 
the CORE sample (bit QSO_EBOSS_CORE) yields the tar¬ 
get density of 50 deg“^ indicated in Tab.[^ 

3.3. Spectroscopic data 

The eBOSS footprint overlaps Stripe 82 over a total of 
120 deg^ delimited by —3° < aj 2 ooo < 45° and —1.25° < 
< 5 j 2000 < 1-25°. The first year of eBOSS observations led to 
the coverage illustrated in Fig.|^ dots indicate the position 
of quasars in plates that have been observed by eBOSS. 
Some outskirt regions (shown in gray in the figure) have 
less than 100 % completeness because of overlapping plates 
that are yet to be observed at the time of this work. In the 
present analysis, we restrict ourselves to the region that has 
100% completeness (in black in the figure). Its total area is 
94.5 deg^. 



R.A. (J2000) 


Fig. 4: Footprint of the first year of eBOSS observations 
used for this study (aspect ratio is not 1:1). Incomplete re¬ 
gions are shown in gray, fully observed ones in black. 


training of the NN was done using a large sample of 13 063 
spectroscopically confirmed quasars with redshifts in 
0.05 < z < 5.0 and magnitudes in the range 18 < g < 23, 
and a star sample consisting in 2 609 objects spectro¬ 
scopically confirmed as stars in the course of the BOSS 
project. The two samples are located in Stripe 82, and 
thus have identical time sampling characteristics as the 
present data. An output unn of the neural network near 
0 designates non-varying objects, as is the case for the 
vast majority of stars, while an output near 1 indicates 
lightcurves exhibiting quasar-like variability. 

A source is selected according to its quasar-like variabil¬ 
ity (unn) and morphology (mdis). A loose morphology cut 
(mdiff < 0 . 1 ) is applied if the variability indicator is high 


In the course of the BOSS survey, several ancillary 
projects have covered parts or all of Stripe 82, either to 
test target selection techniques (Paper Var from quasars 
targeted with BOSS chunk 11 ov er all Stripe 82), or as pilot 


programs for eBOSS and DESI (Ross et al. 2012 Dawson 
et al.|2013 Alam et al.|2015 ). In particular, three programs 
aimed to provide an exnaustive census of quasars to at least 
<7dered ~ 22.5 in Stripe 82. The first one, conducted jointly 
in BOSS (chunk 21, target bit QSO_VAR_FPG) and on the 
Multiple Mirror Telescope, was done in the region delimited 
by 317° < aj 2 ooo < 330°. It provided a total of nearly 1800 
quasars at a mean magnitude (^dered) = 21.1 and led to the 
QLF paper mentioned previously (Paper LF). The other 
two programs were conducted in the region of Stripe 82 de¬ 
limited by 36° < aj 2 ooo < 42° wher e Galactic extinct ion is 
low (c.f., DR12 release of SDSS-III: Alam et al.||2015 ): one 
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program (BOSS chunk 205, target bit QSO_VAR_LF) iden¬ 
tified about 1600 new quasars to guered ~ 22.5, while the 
second (BOSS chunk 218, target bit QSO_DEEP), aimed 
at identifying fainter targets and provided an additional 
363 quasars at (guered) = 22.6. The QSO_VAR_LF and the 
QSO_DEEP programs more than doubled the density of 
known quasars given in Tab. [T] for the rest of Stripe 82. 
The 36° < Q!J 2 ooo < 42° region of Stripe 82 where a deep 
quasar sample is available is used in this work to cross-check 
the completeness-corrected counts computed in Sec. 


4. Quasar number counts 

We here compute the completeness corrections that affect 
our sample, whether they are related to the analysis tech¬ 
nique or to observational constraints. We present raw num¬ 
ber counts derived from the observation of our targets, and 
compute corrected number counts used in Sec. to derive 
a quasar luminosity function. As mentioned above, we only 
give counts within the fully observed zone (black area in 
Fig. 1^, of area 94.5 deg^. 


4.1. Completeness corrections 

We derive the spectroscopy-related completeness correc¬ 
tions from the data themselves, and the selection-related 
corrections from the application of the same selection cuts 
to the control sample of quasars. We describe below the 
different contributions. 


Morphology completeness, emorph(^),ff: 

As explained in Sec. |3.2[ we only select targets among 
the sources that pass the morphology cut mditr < 0.1. 
Therefore, sources that are more extended than allowed 
by this cut are not considered as possible candidates. 
Nevertheless, some low-redshift quasars for which the host 
galaxy is resolved can fail this criterion (c.f.. Fig. [^. We 
compute the correction related to this incompleteness by 
considering the fraction of the quasars in our control sam¬ 
ple that have Wdiff > 0.1, as a function of redshift and 
magnitude. This procedure slightly underestimates the ef¬ 
fect, since the control sample is dominantly consisting of 
objects that were already selected to be unresolved sources. 
The selection, however, was based on single-epoch photom¬ 
etry (and not on coadded images as for the present study), 
making it less sensitive to morphology than the current one. 
The morphology incompleteness 1/emorph is a factor ^1.7 
at redshift z < 0.5, a factor 1.2 for z in 0.5 — 0.8, 1.1 for z 
in 0.8 — 1.0 and is compatible with 1.0 for z > 1. Note that 
we start our measurement of the quasar LF at z = 0.67, 
where the correction is at most of order ^10%. 


Target selection completeness, ese\{g, Wdiff): 

The selection completeness is determined from the fraction 
of quasars in the control sample that pass the color and 
variability criteria of Sec. 3.2 The result is illustrated in 
Fig. [5] as a function of magnitude g and morphology m-diff. 
The efficiency drops for fainter objects where the variability 
signal is not as visible, and for large Wdiff because of the 
stricter variability cut applied to more extended objects. On 


average over the selected sample, the selection completeness 
is 0.86. 



0.9 

0.8 

0.7 

0.6 

0.5 

0.4 


Fig. 5: Selection completeness esei(ff, rudiff) for quasars used 
in this study. The loss at large PSF-model magnitudes, i.e., 
more extended objects, is due to the more stringent cut on 
variability for such sources. The loss at large magnitudes is 
due to increased photometric dispersion in the lightcurves, 
which blurs the variability signal. 


Spectroscopic completeness, espect(5): 

Some spectra did not produce a reliable identification of 
the source, either because the extraction procedure failed 
(yielding flat and useless spectra) or because the spectrum 
had too low a signal-to-noise ratio (S/N) for adequate iden¬ 
tification, whether at the pipeline or at the visual inspection 
level. As explained in Sec.[^ we globally consider such spec¬ 
tra as inconclusive. We make the assumption that the ratio 
of quasars to non-quasars in the identified and inconclusive 
sets are identical. This hypothesis is confirmed by the fact 
that the ensemble structure function parameters (globally 
accounted for by the unn parameter) are similar for the 
two sets. There is a small trend for a higher fraction of 
non-quasars in the inconclusive sample as the magnitude in¬ 
creases beyond g ^ 22, but the effect is small. Furthermore, 
because the selection completeness drops much faster with 
magnitude than the spectroscopic completeness (cf. Fig. 
for a comparative illustration), even a large change in the 
fraction of quasars in the inconclusive set (for instance from 
the current 80% at g ~ 22.5 to 50%) only affects the overall 
completeness at g ~ 22.5 by less than 3%. 

We compute the spectroscopic completeness from the 
fraction of inconclusive spectra as a function of magnitude 
(cf. Fig. [^. This correction is only applied to new quasars 
since previously known ones are, by definition, spectroscop¬ 
ically identified. The completeness correction is 0.93 on av¬ 
erage over the 7900 new quasars, and 1 by definition for 
known quasars. As illustrated in Fig.|^ there is a constant 
~ 1% fraction of inconclusive spectra at bright magnitudes. 
This fraction increases to 8% at g = 22, 16% a.t g = 22.5 
and 30% at g = 23. By comparison, the measured fractions 
of inconclusive spectra are 3%, 8% and 24% at g of 22, 
22.5 and 23, respectively, in the BOSS-I-MMT program of 
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Paper LF. Part of the difference can be explained by the 
fact that we are now more conservative in the identification 
procedure. With identical definitions of inconclusive spec¬ 
tra, the percentages are compatible for g > 22.5, but remain 
slightly larger for the current analysis at g < 22. Inspection 
of the inconclusive spectra at bright magnitudes indicates 
that most of them are at redshift around 0.7 where the au¬ 
tomated pipeline is not yet fully optimized (focus for BOSS 
was on z > 2.2 quasars, and for eBOSS on z > 0.9 quasars, 
thus higher redshifts than where this artefact appears). 



Fig. 6: Fraction of inconclusive spectra (i.e., 1 —espect)^)^) as 
a function of observed g magnitude. The curve is an empir¬ 
ical fit to the data by a hyperbolic tangent. The ~ 1% loss 
at bright magnitudes is an artifact of the current pipeline. 
The increase at faint magnitudes is compatible with previ¬ 
ous estimates. 


Tiling completeness, Ctiiing: 

Some spectroscopic targets cannot be observed either by 
lack of an available fiber on the plate (target density locally 
too high), or because a fiber cannot be placed at that plate 
position, e.g., because of fiber collisions (two fibers cannot 


be located less than 62" apart, c.f. Dawson et al. 20131. 


This loss is random, independent of redshift or magnitude 
and Ctiiing simply indicates the fraction of targets that were 
assigned fibers. It is equal to 0.959 for all new quasars, and 
equal to 1 by definition for already known quasars. 


4.2. Raw number counts 

We identified 7900 new quasars, and selected another 5976 
that had been previously spectroscopically identified. The 
magnitude and redshift distributions of the quasars selected 
by this study are shown in Fig. as the open dark green 
triangles for the new quasars, and as the plain green trian¬ 
gles for the total sample of selected quasars includ ing t he 
already identified ones. The deep sample (cf.. Sec. 3.3) is 


the only one to include quasars beyond gdered = 22.5. 

As is clearly visible from Fig. the newly identified 
quasars have a similar redshift distribution as the previ¬ 
ously identified ones, but extend to fainter magnitudes on 
average, with (t/dered) = 21.5 / 20.4 for the newly / previ¬ 
ously identified quasars, respectively. Tab. lists the raw 
quasar counts for this study, in three bins of observed mag¬ 
nitude g. Our sample is particularly valuable at faint magni¬ 


tudes: at g > 22 , it increases the number of known quasars 
in the footprint by almost a factor of 6 . 


Observed g magnitude 


Sample 

< 21 

21 - 22 

> 22 

Total 

eBOSS 

1206 

4174 

2520 

7900 

Known 

4155 

1389 

432 

5976 


Table 2: Raw number counts for the samples of new 
(‘eBOSS’) and previously identified (‘Known’) quasars for 
the variability-selection of this study, in several g magni¬ 
tude bins. 


4.3. Corrected number counts 


We derive corrected quasar number counts from raw num¬ 
ber counts by accounting for the different sources of in¬ 
completeness detailed above. New quasars are corrected for 
morphology cut, target selection, tiling losses and spectro¬ 
scopic failures, while previously known quasars are only 
corrected for morphology and selection since their iden¬ 
tification does not depend on eBOSS observation. The 
completeness-corrected number of quasars is thus given by 


Nqso 


E 

A^eBOSS 


1 

^sel(,^5 ^diff) ^niorph (-2^ 5 ^) ^tiling Aspect (^) 


E 

Nkt 

^ * K.nown 


1 

^selid^ ^diff) eniorph(.2:, ^) 


(3) 


The total completeness correction for new eBOSS quasars 
has an average of 0.70 and a standard deviation of 0.15, 
and for previously known quasars an average of 0.90 and 
a standard deviation of 0.12. The overall completeness of 
the full sample as a function of magnitude is illustrated in 
Fig. ID 

To validate the computation of the completeness correc¬ 
tions, it is interesting to compare the corrected quasar den¬ 
sities to the measurements done in the deep zone. Despite 
larger counts, it is clearly visible in Fig. (top plot) that 
the deep region still suffers from significant incompleteness 
at 5 dered > 21.5. In Tab. we therefore provide both the 
completeness-corrected densities to a limiting magnitude 
<7dered < 21.5 where the deep sample is expected to be com¬ 
plete, and those to gdered < 22.5 from which we measure 
the quasar LF. At gdered < 21.5, the corrected counts and 
the deep sample show excellent consistency, with less than 
1.5 cr deviations over the entire redshift range (Fig. mid¬ 
dle plot), compatible with Poisson errors. The magnitude 
distributions of the completeness-corrected and deep-zone 
counts (Fig.[^ top plot) are also in excellent agreement to 
5 dered < 21.5. Extending the study to gdered < 22.5 (Tab.[^ 
or Fig. bottom plot), the corrected counts are of order 
10 % larger than the counts in the deep zone, a reasonable 
excess that is easily accounted for by the incompleteness of 
the deep sample at the faint end. 

Tab. also lists the expected number counts from the 
luminosity function computed in Paper LF. The results 
from this new study show a small global increase of 10 
to 20% in number counts over these previous estimates. As 
explained above, the completeness-corrected densities mea¬ 
sured in this work are in agreement with the deep-zone 
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Source 

2 ; < 0.9 

Redshift range 
0.9 < 2 < 2.1 

2 > 2.1 

Total 

Pdered ^ 21.5 
This work 

16 ±1 

64 ±1 

26 ±1 

106 ±1 

Deep zone 

18 ± 1 

60 ±2 

29 ±2 

108 ±3 

Paper LF 

16 

50 

21 

87 

Pdered 22.5 
This work 

23 ±1 

119 ±1 

56 ±1 

198 ±1 

Deep zone 

23 ± 1 

101 ±3 

51 ±2 

175 ±4 

Paper LF 

27 

96 

48 

171 


Table 3: Areal densities (in deg“^) in several redshift bins 
for this work (after completeness correction), for the deep 
36° < a < 42° zone (raw counts) and from the best-fit lu¬ 
minosity function of Paper LF. Uncertainties on raw counts 
are Poisson noise. 


raw counts at the bright end, and show a 10% understand¬ 
able excess at the faint end. Furthermore, the deep sample 
gives a solid lower bound on quasar densities. For these 
reasons, we are confident that the discrepancy between the 
two studies is due to slightly underestimated completeness 
corrections in Paper LF, rather than overestimated correc¬ 
tions in the present work. Given that the total corrections 
in Paper LF ranged from 1.2 at g < 20 to a factor 2 at 
g ~ 22.5, a 10% inaccuracy is not surprising. 


5. Luminosity function in g 

We derive the QLF in a similar manner as described in 
Paper LF. We compute a binned QLF from the corrected 
number counts of Sec. |4.3[ considering our completeness 
limit at gdered = 22.5. We fit two parametric models to our 
binned QLF and compare the number counts each model 
predicts over the range of magnitude and redshift observ¬ 
able by eBOSS. Finally, we use our QLF fits to predict num¬ 
ber counts to fainter magnitudes than achieved by eBOSS, 
as needed for future quasar surveys. 


5.1. Binned luminosity function 

Selection for this survey was performed in the g-band, 
which provides the highest S/N for a vast fraction of the 
data. For each quasar, we compute the absolute magnitude 
normalized to z = 2 by 


Mg{z = 2) = gdered “ dmi^z) - [K{z) - K{z = 2)] 


( 4 ) 


where the distance modulus dM{z) is computed assum¬ 
ing a flat ACDM cosmology with {QA,^M,w ,h) = 


(0.6935,0.3065, —1,0.679) as measured by the Planck 
Collaboration et al. (2015) in the ‘TT-LlowP-l-lensiiig+exF 
configuration, and K{z) is the AT-correction that accounts 
for redshifting of the bandpass of the spectrum. We choose 
to normalize the magnitudes a,t z = 2 because it is close to 
the median redshift of our sample, an d it allows backward 
compatibility with previous studies 


2009 


__ Groom et al. _ 

Ross et al. 2013 Paianque-Delabrouill'e et al.||2013p. As iii 
Paper LF, we use the A-correction as a function of red¬ 


shift that was derived by McGreer et al. |20^) following 
a similar approach as in Richards et al. (2009p. As illus¬ 


trated in Fig . the AT-correction is close to that of| Croorn| 
et al.| ( 2009| for z < 3, but extends to higher redshifts. 


It varies with luminosity, due to the accounting of strong 
quasar emission li nes whose equiv alent widths are a func¬ 
tion of luminosity (Baldwin 1977). This luminosity depen¬ 
dence introduces a spread of ~ 0.25 mag at 2 ^ 2 — 3 where 
the Lyman-a and C-IV lines contribute substantially to the 
flux in g-band. 

We define eight redshift bins, with limits 0.68, 1.06, 1.44, 
1.82, 2.2, 2.6, 3.0, 3.5 and 4.0. The binned LF is computed 
for each r edshift using the mo del-weighted estimator d* sug¬ 
gested by Mi^aji_e^_^ (2001) and used in previous studies 
such as in|Groom et al. (|2009) or in Paper LF. The binned 
LF is given by 




Tyobs 

jS^model 


( 5 ) 


where Mg. and Zi are, respectively, the absolute magnitude 
and the redshift at the center of bin i, is the model 

LF estimated at the center of the bin, is the num¬ 

ber of quasars in the bin with gdered <22.5 estimated from 
the model, and Af°^® is the observed number of quasars 
in the bin. This estimator is free from most of the biases 
unavoidable in the more usual 1/U method devised by ?, 
and improves upon it in several ways. It corrects for vari¬ 
ations of the LF within a bin (particularly critical at the 
steep bright end of the QLF), corrects for incompleteness 
in a bin (particularly critical at the faint end of our LF 
where the bin is incompletely sampled), and allows exact 
errors to be evaluated using Poisson statistics. A dr awback 
of this estim ator is that it is model-dependent, but |Miy^ 
et al.| (|2001|) demonstrated that the uncertainties due to 
the rriodel dependence are practically negligible. We here 
free ourselves of any model-dependence by performing an 
iterative fitting to determine the binned QLF from Eq. 
and our choice of model, until parameter convergence is 
reached. As explained in the following section, we use two 
different models to fit the QLF, which do not produce any 
significant difference to the estimate of the binned LF. 


5.2. QLF model fits 

The QLF is traditional l y fit by a double pow er law of the 
form ( Boyle et al.||2000 Richards et al.|[2ro6 ): 


^Mg,z) = 


$* 


2Q0.4(a+l)(Mg-M;) ^Q0.4(/3-Hl)(Mg-M* 


( 6 ) 


where $ is the quasar comoving space density and M* a 
characteristic, or break, magnitude. The slopes a and /3 
describe the evolution of the LF on either side of the break 
magnitude. In this work, we consider two extensions of this 
simple form described below. 

Our first model is the same as the one we used in 
Paper LF. We consider a pure lu minosity-evolution (PLE) 

where a redshift depen- 


model as in Groom et al. 


(2009) 


dence of the luminosity is introduced through an evolution 
in Mg described by 


Mg{z) = Mg{zp) - 2.5[fci(z - Zp) + k2{z - Zp^] . 


( 7 ) 


We allow the redshift-evolution parameters (fci and /C 2 ) and 
the model slopes (a and jd) to be different on either side 
of a pivot redshift Zp = 2.2. The model is thus described 
by Eqs. and where a, /3, fci and ^2 are defined with 
subscript I for z < Zp and h for z > Zp. This PLE model 
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therefore has ten parameters ($*, M*{zp), ai, Pi, ku, k 2 i, 
oih, Ph, kih and k 2 h) that are free to vary in the fit. 

A n exten sive study of the QLF was performed by |Ross| 


et al. (20131 on 23,300 quasars with i < 2 1.8 and 2.2 < z < 
3.5 from the DR9 release of BOSS data ( Ahn et ^|2012 1, 
complemented by about 5500 quasars over 2.2 < z < 3.5 
from Paper Var and about 1900 quasars over 0.3 < z < 
3.5 from Paper LF. The authors showed that a good fit 
to this large sample of quasars over the full redshift range 
was obtained by using a PLE model for z < 2.2, and a 
model with both luminosity and density evolution (LEDE) 
for 2 > 2.2, where the normalization and break magnitude 
evolve in a log-linear manner, e.g.. 


log[$*(z)] = \og[^*{Zp)] + Cia{z - Zp) + Cib{z - Zpf (8) 


M:{z) = M*{Zp) + C 2 {z-Zp) 


( 9 ) 


with Zp = 2.2 the pivot redshift. This PLE (for z < 2.2) 
-I- LEDE (for z > 2.2) is o ur second model w ith which we 
fit our binned QLF. Unlike Ross et al. (2013), however, we 
impose continuity of the LF at z = Zp by requiring the 
same normalization <i>*(zp) and break magnitude M*{zp) 
for both the PLE and the LEDE forms. We allow for some 
additional flexibility by allowing a redshift dependence of 
the slope, according to 


i{z) = a{zp) -I- C3(z - Zp) 


( 10 ) 


where a{zp) is equal to the value of a used at z < Zp in 
the PLE form. Our PLE-I-LEDE model therefore has ten 
parameters ($*(0), M*(0), a{Zp), P{zp), ki, k 2 , Cia, cu, C 2 
and C3) that are left free to vary in the fit. 

The best-fit parameters are given in Tab. Both mod¬ 
els have ten parameters free to vary in the fit, and start with 
100 data points spread over eight redshift bins. Throughout 
the iterations (a total of ten iterations is needed in both 
cases to reach parameter convergence), the QLF is recom¬ 
puted according to Eq. and points that are corrected 
by more than a factor 2 ^ue to the incompleteness intro¬ 
duced by the ^dered < 22.5 cut) are removed from the fit. 
One to two points per redshift bin are excluded by this 
procedure. The resulting best-fit models are illustr ated in 
Fig. [Tol alo ng with the best-fit model obtained by |Croom 
et al. (2009) for 2: < 2.2 and ex trapolated to high er redshift, 
and the PLE-I-LEDE model of Ross et al. (2013) shown for 
z > 2.2, where it is best constrained by the DR9 data used 
for the fit. In the z <2.2 range, our two models are in excel¬ 
lent agreement with Groom et al. (2009), and are a good fit 
to our binned QLF. At z > 2.2, our two models show simi¬ 
lar trends at the bright end, but start to differ at the faint 
end of the QLF, in particular for the highest two redshift 
bins, which are la cking faint quasars. The agreement with 
Ross et al. (2013) is good over the common redshift range, 
but the fits also kart to deviate at z > 3.5 where data are 
scarce. Although we constrain our models to be continuous 
at z = 2.2, the fit reduced y^s are less than 2. Such values 


were only obtained in Ross et al. 


(2013) when fitting over 


restricted redshift ranges, typically limiting to data below 
or above a redshift of 2.2. Fig. [TT] presents the redshift evo¬ 
lution of the QLF in a series of luminosity bins, including 
both our data and the best-fit PLE-I-LEDE model. The 
‘kink’ at z = 2.2 is due to the change of analytical form 
at this pivot redshift. The model, however, is continuous at 
z = 2.2. 


We provide, in Appendix A, the measured luminos¬ 
ity function values and associated uncertainties, where the 
measurements were corrected using the model-weighted es¬ 
timator of Eq. and considering the best-fit PLE-I-LEDE 
model. 


5.3. Predicted number counts 


Out best-fit luminosity functions can be used to estimate 
the density of quasars as a function of redshift and mag¬ 
nitude. As shown in Tab. the redshift distributions pre¬ 
dicted by our two models agree within 0.3% for counts to 
a limiting magnitude g < 22.5, and within 1.1% to g < 23. 
The largest discrepancy between the models occurs near 
z = 2.3 (cf. Fig. [T^, close to the pivot redshift Zp = 2.2 
where the analytical form of our models changes, causing a 
small discontinuity in the derivative of their redshift evolu¬ 
tion. At z ^ 2.3, the predictions from the two models differ 
by about 10%. 

The two models are also in goo d agreement with the 
corrected number counts of Sec. 4.3 (Tab. [^, and indicate 
a ~ 15% increase over previous estimates based on the QLF 
fit of Paper LF. Both the PLE and the PLE-I-LEDE mod¬ 
els ht the z > 2 data well, but do not provide as good a fit 
to the z < 0.6 range. A possible explanation for this dis¬ 
crepancy is a loss of low-z quasars in our selection, possibly 
due to the morphology cut for which the completeness ef¬ 
ficiency Emorph does not sufficiently account. Another pos¬ 
sibility is the misidentification of low-z quasars with the 
current pip eline (cf. paragraph on spectroscopic complete¬ 
ness in Sec. 4.1). This small discrepancy, however, is of lit¬ 
tle relevance for large-scale spectroscopic surveys that are 
mostly focusing on quasars at z > 1 where they are abun¬ 
dant enough to probe matter clustering. 


Model 

Limit 

Redshift range 
<0.9 0.9-2.1 >2.1 

Total 

PLE 

g < 22.5 

32 

111 

53 

196 

PLE-hLEDE 

g < 22.5 

31 

108 

55 

195 

PLE 

5 < 23 

43 

151 

75 

269 

PLE-hLEDE 

5 < 23 

41 

145 

77 

263 

PLE 

r < 23 

49 

171 

91 

311 

PLE-^LEDE 

r <23 

46 

162 

88 

296 


Table 5: Quasar densities (in deg“^) predicted from the two 
phenomenological models of this study. The PLE and the 
PLE-I-LEDE models give very similar density estimates. 


Using the average g — r vs. redshift dependence that we 
measure for all DR12Q BOSS quasars, we can also calculate 
the number of quasars as a function of r-band magnitude. 
We provide these estimates for our two models as a function 
of g in Tab. and as a function of r in Tab. [Tj for an 
hypothetical survey covering 10,000 deg^. 

We can apply our best-ht QLF to future surveys, such as 
the third-generation large-scale spectroscopic survey DESI 
that aims to observe quasars to a limiting magnitude g = 
23, or possibly r = 23 (in order to recover quasars at z > 3.6 
that are 5-band dropouts due to the absorption of their flux 
by Lyman-a absorbers along the line of sight). Expected 
quasar densities at DESI depth for our two models are given 
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Model 

Redshift 

range 



Parameters 



X^lv 

PLE 

0.68-4.0 

0.68 - 2.2 
2.2-4.0 

Mg (Zp) 

—26.71±o.i5 

log($‘) 

—6.01±o.o7 

a 

—4.31±0.26 
—3.04±o.i2 

p 

— 1.54±o.04 

— 1.38±o.o7 

ki 

—0.08±o.o8 
—0.25±o.o9 

k2 

—0.40±o.05 
—0.05±o.o6 

135/76 

PLE 
+ LEDE 

0.68-4.0 

0.68 - 2.2 

m;(o) 

—22.25±o.49 
ki 

1.59±0.28 

log[4>*(0)] 

—5.93±o.o9 
k2 

—0.36±o.o9 

a 

—3.89±o.23 

P 

—1.47±o.o6 



146/77 


2.2-4.0 

Cla 

—0.46±o.io 

Cu 

—0.06±o.io 

C2 

—0.14±o.i7 

C3 

0.32±0.23 





Table 4: Values of the parameters (and redshift ra^e over which they apply) for the best-fit PLE and PLE-I-LEDE 
models of quasar luminosity functions (e.g., Eqs. [6[J^. The slope a reproduces the bright end part of the QLF, and /3 
the faint end. 



-29 -28 -27 -26 -25 -24 -23 -22 -21 I9 -28 -27 -26 -25 -24 -23 -22 -21 I9 -28 -27 -26 -25 -24 -23 -22 -21 !9 -28 -27 -26 -25 -24 -23 -22 -21 

Mg(z=2) Mg(z=2) Mg(z=2) Mg(z=2) 


Fig. 10: Quasar luminosity function measurements (black circles). The best-fit models of this work are shown as the red 


(respectively blue) curves for the PLE-I-LEDE (resp. PLE) models. The green dot-dashed curve is the LF of Groom et al. 
(2009). The plain cyan curve is the best fit LEDE model of|Ross et ah' (2013) at z > 2.2. The orange dotted and dashed 
curves are the best fits to COSMOS data (Masters et al. 2012 [ ) at z ~ 3.2 (shown in the last two redshift bins) and 
z ~ 4.0, respectively. The magenta dashed curve (almost exactly overlapping th e orange dotted cu rve at the faint end in 
the 3.0 < z < 3.5 redshift bin) is measured at z ~ 3.2 from SWIRE and SDSS (Siana et al.||2008). 


in the last two sections of Tab. [5] The 0.9 < z < 2.1 redshift 
range is where quasars are currently used as direct tracers 
of dark matter, and the z > 2.1 regime is where quasars 
are used to probe dark matter though the Lyman-a forest. 
The r < 23 limit produces slightly higher estimates than 


at g < 23, because quasars have positive g — r values at all 
redshift {g — r in 0.1 - 0.5 for z < 3.5, and increasing at 
higher redshift). Despite the differences between the PLE 
and the PLE-I-LEDE models that are visible in Fig.[^or in 
the predicted counts of Tabs. andfor faint magnitudes in 
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g\z 

0.5 

1.5 

PLE model 

2.5 3.5 

4.5 

5.5 

Total 

15.75 

26 

2 

10 

0 

0 

0 

38 

16.25 

58 

8 

33 

0 

0 

0 

100 

16.75 

143 

39 

87 

3 

0 

0 

272 

17.25 

386 

180 

231 

10 

0 

0 

806 

17.75 

1134 

813 

626 

26 

0 

0 

2599 

18.25 

3396 

3468 

1727 

68 

0 

0 

8660 

18.75 

8930 

12525 

4780 

179 

0 

0 

26413 

19.25 

17963 

32977 

12319 

462 

1 

0 

63722 

19.75 

28172 

61064 

26357 

1170 

3 

0 

116766 

20.25 

38732 

90121 

45742 

2827 

7 

0 

177430 

20.75 

50756 

121102 

68898 

6299 

18 

0 

247074 

21.25 

65550 

157999 

95351 

12540 

47 

0 

331487 

21.75 

84300 

204485 

125298 

22098 

120 

0 

436302 

22.25 

108276 

264163 

159894 

35046 

301 

1 

567680 

22.75 

138990 

341219 

201144 

51351 

737 

3 

733442 

23.25 

178321 

440913 

251597 

70978 

1722 

7 

943538 

23.75 

228612 

570015 

314273 

93793 

3739 

17 

1210449 

24.25 

292750 

737314 

392821 

119723 

7357 

42 

1550006 

24.75 

374172 

954240 

491772 

149132 

13010 

106 

1982433 

Total 

1620667 

3992648 

2192958 

565706 

27062 

175 

8399216 



PLE-I-LEDE model 




9\z 

0.5 

1.5 

2.5 

3.5 

4.5 

5.5 

Total 

15.75 

35 

5 

1 

0 

0 

0 

42 

16.25 

82 

20 

5 

0 

0 

0 

108 

16.75 

200 

77 

20 

1 

0 

0 

298 

17.25 

519 

291 

76 

3 

0 

0 

889 

17.75 

1402 

1088 

287 

10 

0 

0 

2787 

18.25 

3758 

3887 

1065 

35 

0 

0 

8745 

18.75 

9092 

12345 

3737 

117 

0 

0 

25291 

19.25 

17972 

31289 

11415 

382 

1 

0 

61059 

19.75 

28629 

59594 

27445 

1163 

4 

0 

116836 

20.25 

39554 

90405 

50511 

3077 

13 

0 

183560 

20.75 

51163 

121609 

76197 

6648 

36 

1 

255653 

21.25 

64559 

155923 

103193 

11780 

97 

2 

335554 

21.75 

80771 

196666 

133441 

18077 

242 

4 

429201 

22.25 

100758 

246795 

169698 

25349 

543 

9 

543153 

22.75 

125546 

309293 

214695 

33716 

1063 

21 

684334 

23.25 

156321 

387578 

271321 

43594 

1820 

46 

860680 

23.75 

194486 

485819 

342970 

55625 

2799 

93 

1081792 

24.25 

241694 

609219 

433845 

70598 

3989 

173 

1359518 

24.75 

299797 

764319 

549267 

89454 

5393 

291 

1708521 

Total 

1416338 

3476223 

2389189 

359629 

16003 

640 

7658021 


Table 6: Predicted differential quasar counts over 15.5 < g < 25 and 0 < z < 6 for a survey covering 10,000 deg^, based on 
our best-fit PLE or PLE-I-LEDE luminosity function model. Bins are centered on the indicated magnitude and redshift 
values. The ranges in each bin are Ag = 0.5 and Az = 1. 


particular, the density of quasars predicted by either model 
over the redshift and magnitude ranges of next-generation 
surveys are in excellent agreement, at the 1.1% level for 
g <23 and the 2.5% level for r < 23 in total quasar counts. 

6. Conclusions 

The extended Baryon Oscillation Spectroscopic Survey 
(eBOSS), part of the fourth iteration of the Sloan Digital 
Sky Survey, has an extensive spectroscopic quasar program 
that combines several selection techniques. Algorithms us¬ 
ing the variability of quasar luminosity with time have been 
shown to be highly efficient to obtain large and complete 
samples of quasars. 

Here we present the use in eBOSS of time-domain vari¬ 
ability, and focus on a specific program in Stripe 82 that 


led to a sample of 13,876 quasars to ^dered = 22.5 over a 
94.5 deg^ region, 1.5 times denser than expected to be ob¬ 
tained over the rest of the eBOSS footprint. This variabil¬ 
ity program provides a homogeneous and highly complete 
sample of quasars that is denser and of greater depth than 
samples that were used in previous studies dedicated to 


QLF measurements such as 

Groom et al. 

to 

o 

o 

Palanque- 

Delabrouille et al. (2013); 

Ross et al.| (|2013|). Using the 


data themselves, plus an external control sample of quasars 
selected with an independent (color-based) technique, we 
compute completeness corrections to account for quasar 
losses in the selection procedure, in the tiling, or in the 
spectroscopic identification. 

We use this sample to measure the QLE in eight red¬ 
shift bins from 0.68 to 4.0, and over magnitudes ranging 
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r\z 

0.5 

1.5 

PLE model 

2.5 3.5 

4.5 

5.5 

Total 

15.75 

54 

5 

14 

0 

0 

0 

73 

16.25 

127 

23 

41 

3 

0 

0 

194 

16.75 

323 

104 

108 

9 

0 

0 

544 

17.25 

860 

469 

290 

25 

1 

0 

1644 

17.75 

2316 

2009 

793 

64 

2 

0 

5185 

18.25 

5854 

7439 

2204 

167 

5 

0 

15669 

18.75 

12619 

21145 

5994 

432 

12 

1 

40203 

19.25 

22090 

44412 

14675 

1100 

31 

2 

82311 

19.75 

32538 

72781 

29834 

2700 

80 

5 

137938 

20.25 

43787 

102680 

50203 

6208 

201 

12 

203090 

20.75 

57009 

136093 

74188 

12952 

502 

30 

280774 

21.25 

73477 

176857 

101326 

24039 

1217 

75 

376990 

21.75 

94414 

228620 

132051 

39635 

2811 

189 

497720 

22.25 

121186 

295216 

167766 

58919 

6038 

470 

649595 

22.75 

155438 

381240 

210616 

80786 

11825 

1144 

841048 

23.25 

199196 

492527 

263236 

104700 

20944 

2660 

1083262 

23.75 

254945 

636610 

328743 

131042 

33694 

5769 

1390804 

Total 

1076235 

2598227 

1382081 

462783 

77360 

10356 

5607043 



PLE-I-LEDE model 




r\z 

0.5 

1.5 

2.5 

3.5 

4.5 

5.5 

Total 

15.75 

74 

12 

2 

0 

0 

0 

88 

16.25 

173 

46 

8 

1 

0 

0 

227 

16.75 

425 

173 

29 

3 

0 

0 

630 

17.25 

1063 

645 

109 

11 

1 

0 

1830 

17.75 

2637 

2318 

407 

37 

3 

1 

5404 

18.25 

6171 

7534 

1476 

119 

10 

3 

15314 

18.75 

12760 

20293 

4934 

380 

27 

7 

38401 

19.25 

22256 

42853 

14012 

1146 

71 

15 

80354 

19.75 

33079 

72019 

31432 

3068 

182 

33 

139813 

20.25 

44377 

102993 

55189 

6846 

432 

69 

209907 

20.75 

56787 

135480 

81153 

12487 

915 

133 

286955 

21.25 

71385 

172381 

108661 

19287 

1689 

232 

373634 

21.75 

89178 

216865 

139912 

26781 

2731 

367 

475834 

22.25 

111154 

271898 

177682 

35182 

3976 

534 

600426 

22.75 

138390 

340627 

224720 

45094 

5393 

732 

754956 

23.25 

172133 

426753 

283998 

57240 

7015 

962 

948101 

23.75 

213829 

534834 

359048 

72418 

8921 

1234 

1190284 

Total 

975871 

2347722 

1482773 

280101 

31368 

4322 

5122156 


Table 7: Predicted differential quasar counts over 15.5 < r < 24 and 0 < z < 6 for a survey covering 10,000 deg^, based on 
our best-fit PLE or PLE-I-LEDE luminosity function model. Bins are centered on the indicated magnitude and redshift 
values. The ranges in each bin are Ag = 0.5 and Az = 1. 


from Mg{z = 2) = —26.60 to —21.80 at low redshift, and 
from Mg{z = 2) = —29.00 to —26.20 at high redshift. The 
data indicate a break at pivot redshift Zp = 2.2, with a rise 
in luminosity followed by a steep decline as the redshift 
increases on either side of Zp. The data are well fit by a 
double power-law model. We compare two models that we 
constrain to be continuous at z^: a, quadratic PLE model as 


Palanque-Delabrouille et al. (20131, with bright-end and 


faint-end slopes allowed to be dif termt on either sid e of Zp, 
and a PLE -|- LEDE model as in Ross et al. (2013]), where 
a simple linear PLE is used for z < 2.2 and a LEDE with 
quadratic magnitude evolution is used at z > 2.2. These 
models both have ten parameters free to vary, and they fit 
the measured binned QL E equally well. Our m odels are in 
excellent agreement wit h Groom et al. (2009) at z < 2.2, 
and with {Masters et al. (2012) at z > 2.2. Our two models 
start to deviate from one another in the highest two redshift 
bins (z > 3.0), although they both are in reasonable agree- 


ment with other measurements (Siana et al.||2008 

Masters 

et al.||2012 

Ross et al.||2013 

)• 


We use our models to predict densities of quasars for fu¬ 
ture quasar spectroscopic surveys. We predict 266±3 deg“^ 
quasars to g < 23, and 304 ± 7 deg” quasars to r < 23, 
where the estimate is the mean of the two model estimates 
and the uncertainty is the difference between the estimate 
from either model and the mean. 


Appendix 

Tab. provides the binned luminosity function measured 
in this work using the model-weighted estimator described 
in Sec. |5.2[ and plotted in Fig. The table lists the value 
of log$ in eight intervals spanning redshifts from z = 0.68 
to z = 4.00, and for AMg = 0.40 magnitude bins from 
Mg = —29.00 to Mg = —20.60. We also give the number 
of quasars {Nq) contributing to the LF in each bin, and 
the uncertainty (Alog$). Bins with quasars but no cor- 
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responding value of the binned QLF are data points that 
were removed in the iterative fitting procedure due to large 
correction factors. 

Acknowledgements. Funding for the Sloan Digital Sky Survey IV 
has been provided by the Alfred P. Sloan Foundation, the U.S. 
Department of Energy Office of Science, and the Participating 
Institutions. SDSS-IV acknowledges support and resources from the 
Center for High-Performance Computing at the University of Utah. 
The SDSS web site is www.sdss.org. 

SDSS-IV is managed by the Astrophysical Research Consortium 
for the Participating Institutions of the SDSS Collaboration includ¬ 
ing the Brazilian Participation Group, the Carnegie Institution for 
Science, Carnegie Mellon University, the Chilean Participation Group, 
the French Participation Group, Harvard-Smithsonian Center for 
Astrophysics, Institute de Astroffsica de Canarias, The Johns Hopkins 
University, Kavli Institute for the Physics and Mathematics of the 
Universe (IPMU) / University of Tokyo, Lawrence Berkeley National 
Laboratory, Leibniz Institut fiir Astrophysik Potsdam (AIP), Max- 
Planck-Institut fiir Astronomie (MPIA Heidelberg), Max-Planck- 
Institut fiir Astrophysik (MPA Garching), Max-Planck-Institut fiir 
Extraterrestrische Physik (MPE), National Astronomical Observatory 
of China, New Mexico State University, New York University, 
University of Notre Dame, Observatario Nacional / MCTI, The 
Ohio State University, Pennsylvania State University, Shanghai 
Astronomical Observatory, United Kingdom Participation Group, 
Universidad Nacional Autonoma de Mexico, University of Arizona, 
University of Colorado Boulder, University of Portsmouth, University 
of Utah, University of Virginia, University of Washington, University 
of Wisconsin, Vanderbilt University, and Yale University. 

We acknowledge D. Vilanova for his help in the selection of the quasar 
targets with the variability algorithm. 


References 

Ahn C.P., Alexandroff R., Allende Prieto C., et ah, Dec. 2012, The 
Astrophysical Journal Supplement Series, 203, 21 
Alam S., Albareti F.D., Allende Prieto C., et al., Jan. 2015, ArXiv 
e-prints 

Annis J., Soares-Santos M., Strauss M.A., et al., Oct. 2014, ApJ, 794, 
120 

Baldwin J.A., Jun. 1977, ApJ, 214, 679 

Bolton A.S., Schlegel D.J., Aubourg E., et ah, Nov. 2012, The 
Astronomical Journal, 144, 144 

Boyle B.J., Shanks T., Peterson B.A., Dec. 1988, MNRAS, 235, 935 
Boyle B.J., Shanks T., Groom S.M., et al., Oct. 2000, MNRAS, 317, 
1014 

Busca N.G., Delubac T., Rich J., et al., Apr. 2013, Astronomy 
Astrophysics, 552, A96, URL http://arxiv.org/abs/1211. 
2616http://WWW.aanda.org/10.1051/0004-6361/201220724 
Groom S.M., Shanks T., Boyle B.J., et al., Aug. 2001, MNRAS, 325, 
483 

Groom S.M., Smith R.J., Boyle B.J., et al., Apr. 2004, MNRAS, 349, 
1397 

Groom S.M., Richards G.T., Shanks T., et ah, Jan. 2009, MNRAS, 
392, 19 

Dawson K., et ah, in preparation, 2015, to be submitted to ApJ 
Dawson K.S., Schlegel D.J., Ahn C.P., et ah, Jan. 2013, The 
Astronomical Journal, 145, 10 

Delubac T., Bautista J.E., Busca N.G., et ah, Feb. 2015, A&A, 574, 
A59 

Doi M., Tanaka M., Fukugita M., et ah, Apr. 2010, AJ, 139, 1628 
Eisenstein D.J., Zehavi L, Hogg D.W., et ah, Nov. 2005, ApJ, 633, 
560 

Eisenstein D.J., Weinberg D.H., Agol E., et ah, Sep. 2011, The 
Astronomical Journal, 142, 72 
Fan X., May 1999, AJ, 117, 2528 

Fan X., Strauss M.A., Schneider D.P., et ah, Jan. 2001, AJ, 121, 54 
Fukugita M., Ichikawa T., Gunn J.E., et ah, Apr. 1996, AJ, 111, 1748 
Gunn J.E., Siegmund W.A., Mannery E.J., et ah, Apr. 2006, The 
Astrophysical Journal, 131, 2332 
Hasinger G., Miyaji T., Schmidt M., Oct. 2005, A&A, 441, 417 
Hopkins P.F., Richards G.T., Hernquist L., Jan. 2007, ApJ, 654, 731 
Koo D.C., Kron R.G., Feb. 1988, ApJ, 325, 92 

Law N.M., Kulkarni S.R., Dekany R.G., et ah, Dec. 2009, PASP, 121, 
1395 

Masters D., Capak P., Salvato M., et ah, Aug. 2012, ApJ, 755, 169 


McGreer I.D., Jiang L., Fan X., et ah. May 2013, ApJ, 768, 105 
Miyaji T., Hasinger G., Schmidt M., Apr. 2001, A&;A, 369, 49 
Myers A., et ah, in preparation, 2015, to be submitted to ApJ 
Palanque-Delabrouille N., Yeche C., Myers A.D., et ah, Jun. 2011, 
A&A, 530, A122 

Palanque-Delabrouille N., Magneville G., Yeche C., et ah, Feb. 2013, 
Astronomy & Astrophysics, 551, A29 
Paris L, Petitjean P., Aubourg E., et ah, Nov. 2012, Astronomy Sz 
Astrophysics, 548, A66 

Paris I., Petitjean P., Aubourg E., et ah. Mar. 2014, A&A, 563, A54 
Planck Collaboration, Ade P.A.R., Aghanim N., et ah, Feb. 2015, 
ArXiv e-prints 

Rau A., Kulkarni S.R., Law N.M., et ah, Dec. 2009, PASP, 121, 1334 
Richards G.T., Fan X., Newberg H.J., et ah, Jun. 2002, AJ, 123, 2945 
Richards G.T., Groom S.M., Anderson S.F., et ah, Juh 2005, MNRAS, 
360, 839 

Richards G.T., Strauss M.A., Fan X., et ah, Jun. 2006, AJ, 131, 2766 
Richards G.T., Myers A.D., Gray A.G., et ah, Jan. 2009, ApJS, 180, 
67 

Ross N.P., Myers A.D., Sheldon E.S., et ah. Mar. 2012, The 
Astrophysical Journal Supph, 199, 3 
Ross N.P., McGreer I.D., White M., et ah, Aug. 2013, ApJ, 773, 14 
Schlegel D., Abdalla F., Abraham T., et ah, Jun. 2011 
Schlegel D.J., Finkbeiner D.P., Davis M., Jun. 1998, ApJ, 500, 525 
Schmidt K.B., Marshall P.J., Rix H.W., et ah, May 2010, ApJ, 714, 
1194 

Schmidt M., Green R.F., Jun. 1983, ApJ, 269, 352 
Schneider D.P., Richards G.T., Hall P.B., et ah, Jun. 2010, AJ, 139, 
2360 

Siana B., Polletta M.d.C., Smith H.E., et ah. Mar. 2008, ApJ, 675, 49 
Slosar A., Irsic V., Kirkby D., et ah, Apr. 2013, Journal of Cosmology 
and Astroparticle Physics, 2013, 026 
Smee S.A., Gunn J.E., Uomoto A., et ah, Aug. 2013, The 
Astrophysical Journal, 146, 32 

Stoughton C., Lupton R.H., Bernard! M., et ah, Jan. 2002, AJ, 123, 
485 

Tinker J., SDSS-IV Collaboration, Jan. 2015, In: American 
Astronomical Society Meeting Abstracts, voh 225 of American 
Astronomical Society Meeting Abstracts, #125.06 
York D.G., Adelman J., Anderson, Jr. J.E., et ah, Sep. 2000, The 
Astronomical Journal, 120, 1579 


13 


N. Palanque-Delabrouille et al.: Quasar Luminosity Function 



(0 

(A 

o 

c 

0) 

Q. 

E 

o 

o 


1 



0.6 - 


0.4 


0.2 - 


Morphology 

Morphology x Selection 

Morphology x Selection x Tiling x Spectro 


18 18.5 


I I I 

19 


I I I I I I I I I 

19.5 20 


L 


I I I 

20.5 21 21.5 22 22.5 

Observed magnitude g 




Fig. 8: Magnitude dependence of the contributions to the 
total completeness correction, averaged over all variability- 
selected quasars of the eBOSS and known samples. 



Fig. 9: The g-band K-correction as a function of redshift. 
The spread illustrates the luminosity variation of the cor¬ 
rection for the quasars in our sample. The sudden change 
at z > 3 is due to the suppression of the flux in g-band by 
the intervening Lyman-a absorbers. 


Fig. 7: Extinction-corrected magnitude (top) and redshift 
(middle and bottom) distributions. At faint magnitudes, 
most variability-selected sample (green triangles) comes 
from the newly-identified quasars (open dark green trian¬ 
gles). The deep sa mple (orange circles) from the 36° < a < 
42° zone (cf. Sec. reproduces well the corrected counts 
(plain red circles) to g = 21.5, validating the computation 
of the completeness corrections to this magnitude limit. The 
deep sample is the only one to extend beyond g = 22.5. 
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Fig. 11: . Comparison of our best-fit PLE-I-LEDE model 
(solid curves) to the eBOSS Stripe 82 QLF data (points). 
From bottom to top, the curves correspond to magnitudes 
Mg{z = 2) = -28.4, -28, -27.6, -27.2, -26.8, -26.4, 
-26, -25.6, -25.2, -24.8, -24.2, -24 and -23.6. A simple 
PLE model breaks down at z = 2.2. Here, we switch to an 
LEDE model at z > 2.2. The mixed PLE-I-LEDE model 
reproduces the QLF well. 



Fig. 12: Projected counts for gc < 22.5 for the PLE (blue) 
and the PLE+LEDE (red) luminosity function models used 
to fit the binned QLF. The dashed curve is the best-fit 
model of Paper LF, shown here for comparison. Our new 
fits indicate a ~ 10% increase in total quasar counts. 
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2.20 < 2 : 

< 2.60 

2.60 < 2 : 

< 3.00 

3.00 < 2 

< 3.50 

3.50 < 2 

< 4.00 

(bin center) 

Nq 

log 4? 

A log 4? 

Nq 

log 4? 

A log 4? 

Nq 

log 4? 

A log 4? 

Nq 

log 4? 

A log 4? 

-28.80 

0 

- 

- 

0 

- 

- 

0 

- 

- 

1 

-8.34 

0.44 

-28.40 

2 

-7.98 

0.31 

7 

-7.42 

0.17 

3 

-7.89 

0.26 

1 

-8.35 

0.44 

-28.00 

6 

-7.49 

0.18 

15 

-7.10 

0.12 

11 

-7.33 

0.14 

1 

-8.36 

0.44 

-27.60 

25 

-6.88 

0.10 

20 

-6.99 

0.11 

11 

-7.33 

0.14 

0 

- 

- 

-27.20 

41 

-6.67 

0.08 

31 

-6.80 

0.09 

13 

-7.26 

0.13 

7 

-7.49 

0.17 

-26.80 

79 

-6.39 

0.07 

52 

-6.57 

0.08 

41 

-6.75 

0.08 

13 

-7.17 

0.13 

-26.40 

125 

-6.18 

0.06 

91 

-6.32 

0.07 

61 

-6.55 

0.07 

25 

-6.71 

0.10 

-26.00 

132 

-6.15 

0.06 

124 

-6.16 

0.06 

85 

-6.37 

0.07 

17 

- 

- 

-25.60 

192 

-5.96 

0.06 

151 

-6.05 

0.06 

105 

-6.19 

0.06 

6 

- 

- 

-25.20 

267 

-5.79 

0.06 

177 

-5.93 

0.06 

69 

-6.23 

0.07 

2 

- 

- 

-24.80 

337 

-5.66 

0.05 

169 

-5.88 

0.06 

53 

-6.10 

0.08 

1 

- 

- 

-24.40 

282 

-5.67 

0.05 

140 

-5.81 

0.06 

6 

- 

- 

0 

- 

- 

-24.00 

230 

-5.62 

0.06 

20 

- 

- 

0 

- 

- 

0 

- 

- 

-23.60 

35 

- 

- 

0 

- 

- 

0 

- 

- 

0 

- 

- 

-23.20 

0 

- 

- 

0 

- 

- 

0 

- 

- 

0 

- 

- 

-22.80 

0 

- 

- 

0 

- 

- 

0 

- 

- 

0 

- 

- 

-22.40 

0 

- 

- 

0 

- 

- 

0 

- 

- 

0 

- 

- 

-22.00 

0 

- 

- 

0 

- 

- 

0 

- 

- 

0 

- 

- 

-21.60 

0 

- 

- 

0 

- 

- 

0 

- 

- 

0 

- 

- 

-21.20 

0 

- 

- 

0 

- 

- 

0 

- 

- 

0 

- 

- 

-20.80 

0 

- 

- 

0 

- 

- 

0 

- 

- 

0 

- 

- 


Table 8 : Bi nned quasar luminosity function. The 13,876 quasars were selected with the variability algorithm presented 
in Sec. 3.2 at gdered < 22.5. All quasars lie in a 94.5 deg^ region of Stripe 82 fully observed by eBOSS during the first 


year of the survey. The corresponding QLF points are shown in Fig. 10 
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